arXiv:1504.01397v2 [astro-ph.CO] 29 Dec 2015 


Draft version December 31, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


IDCS J1426.5+3508: THE MOST MASSIVE GALAXY CLUSTER AT Z > 1.5 

Mark Brodwin\ Michael McDonald^, Anthony H. Gonzalez^, S. A. Stanford'^’®, Peter R. Eisenhardt®, 

Daniel Stern® & Gregory R. Zeimann’’ 

Draft version December 31, 2015 

ABSTRACT 

We present a deep (100 ks) Chandra observation of IDCS J1426.5+3508, a spectroscopically con¬ 
firmed, infrared-selected galaxy cluster at z = 1.75. This cluster is the most massive galaxy cluster 
currently known at z > 1.5, based on existing Sunyaev-Zeldovich (SZ) and gravitational leasing de¬ 
tections. We confirm this high mass via a variety of X-ray scaling relations, including T^-M, /g-M, 
Yx-M, and Lx-M, finding a tight distribution of masses from these different methods, spanning 
M 500 = 2.3-3.3 X 10^"^ Mq, with the low-scatter Yx-based mass M^oo,Yx = ^ ^0- 

IDCS J1426.5+3508 is currently the only cluster at z > 1.5 for which X-ray, SZ and gravitational 
lensing mass estimates exist, and these are in remarkably good agreement. We find a relatively tight 
distribution of the gas-to-total mass ratio, employing total masses from all of the aforementioned in¬ 
dicators, with values ranging from fgas ,500 = 0.087-0.12. We do not detect metals in the intracluster 
medium (ICM) of this system, placing a 2cr upper limit of Z(r < i? 5 oo) < 0.18 Zq. This upper limit on 
the metallicity suggests that this system may still be in the process of enriching its ICM. The cluster 
has a dense, low-entropy core, offset by ^30 kpc from the X-ray centroid, which makes it one of the 
few “cool core” clusters discovered at z > 1, and the first known cool core cluster at z > 1.2. The 
offset of this core from the large-scale centroid suggests that this cluster has had a relatively recent 
(< 500 Myr) merger/interaction with another massive system. 

Subject headings: galaxies: clusters: individual (IDCS J1426.5+3508) — galaxies: clusters: intraclus¬ 
ter medium — galaxies: high-redshift — large scale structure of universe — X-rays: 
galaxies: clusters 


1. INTRODUCTION 

In recent years the study of galaxy clusters has mean¬ 
ingfully entered the z > 1 regime, wi th several surveys 


2011 

Mehrtens et al. 2012 

), Sunyaev-Zel’dovich (SZ, 

Bleem et ai.||2015t IHasselhe 

Id et al.||2013| IPlanck Col- 

laboration et al.|2015|), infrared (IR, Eisenhardt et al. 

2DD5[ 

Aluzzm et al. 1120091 IPapovich et al. 1120101 IRettura 

et al.| 

20141 and radio (IWylezalek et al. 2013 Calametz 

et al.|2U13[ Blanton et ai.|2U14| R. Raterno-Mahler et al. 


tende d the reach of cluster co smology (e.g., Benson et al. 
20131, scaling relations (e.g., Andersson et al.|[2011) and 


galaxy for mation and evolutioh m the richest environ¬ 


ment s (e.g. , Tran et aH2010 Hilton et al.|2010' Brodwin 


et al.||2bl3 ) to z ~ 1.5. 


It is crucial to identify the earliest massive progenitors 
of these 1 ^ 2 ^ 1-5 cluster samples, and the present- 
day massive clusters into which they evolve, in order to 
quantify the build-up of the intracluster medium (ICM) 
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and the establishment of self-similarity. In particular, 
the scaling relations between different cluster mass ob¬ 
servables — e.g., ICM temperature, SZ signal, weak lens¬ 
ing — are calibrated to better than ^ 20% at z < 0.5 
.g., Kravtsov et al.|2006 Vikhlinin et al.|2009b Ander- 


sson et al.|2011 ), but are currently poorly constrained at 


high redshitts (z > 1) where clusters provide the greatest 
leverage as probes of the growth of structure. 

Very massive, high redshift galaxy clusters also provide 
a natural testing ground to confirm or refute the predic¬ 
tions stemming from recent galaxy evolution studies in 
clusters at 1 < z < 1.5. While little star formation ac¬ 
tivity is seen in the most massive (M ^ 10^® Mpi) South 


Pole Telescope 

SPT) clusters in this regime (e.g., 

Brod-1 

win et al. 

2010 

Foley et al. 

201 1| Stalder et al. 



the range M ~ (1 — 4) x 10^"^ Mq. Indeed, Brodwin 


et al.| (2013) reported the discovery of a major epoch of 
star format ion at z > 1.4 for IRAC Shallow Cluster Sur¬ 
vey (ISCS, jEisenhardt et al.||2008|) clusters in this mass 
range. They attributed it in part to a high me rging rate 
of gas-rich cluster members, as suggested by |Mancone| 
eT^ ( [Mot . A consequence of this model is that the 
epoch at which the merging (and hence merger-induced 
star formation) ceases should be a function of cluster 
mass. While the SPT clusters at z < 1.3 are generally 
too massive to permit efficient merging, their lower-mass 
progenitors at higher redshifts should have the high star 
formation rates seen in the ISCS. Indeed, in the most 
distant SPT cluster, SPT-CL J2040-4451, with a mass 
of Msoo.sz 3.2 X 10^"* Mq at z = 1.48, residual star 
formation activity still persists after its earlier bursting 
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phase (Bayliss et al.|[2M4) 


To extend the inlrared cluster search to z > 1.5, we 
repeated the methodology of the IS CS using the Spitze r 
Deep, Wide-Field Survey (SDWFS, [Ashby et al. 2009), 


which quadrupled the Spitzer/I RAC exposure time ove r 
the 9 deg^ IRAC Shallow Survey (Eisenhardt et al.|2004). 


(IDCS, 

Stanford et al. 

2012 

Brodwin et al. 

2012 

Gon- 

zalez et al.||2012D, has identihed two ot the most distant 

clusters to date: IDCS J 1426.5+3508 at z = 

1.75 (Stan- 


1 Zeimann et ai.|2012 ). The latter appears to be a moder- 
ate mass cluster still in the process of formation, whereas 
IDCS J1426.5+3508 at z = 1.75, the subject of this pa- 
p er, is a very ma s sive c luster. 


Stanford et al. 


(2012) reported a detection in only 8.3 
ks of archival Uhandra X-ray Observatory data, result¬ 
ing in an Lx-based mass estimate of M 500 , Lx = (3-3 ± 
1.0) X 10^"^ M(^ IDCS J1426.5+3508 was also observed 
with the Sunyaev-Zeldovich Array (SZA), a subarray 
of the Combined Array for Research in Millimeter-wave 
Astronomy (CARMA). An SZ-based mass of Msoo.sz 
= (2.6 ± 0.7) X 10^^ Min was measured from the strong 


background flares before applying the latest calibration 
corrections using CIAO v4.6 and CALDB v4.6.1.1. 

3. X-RAY PROPERTIES OE IDCS J1426.5+3508 

3.1. Images, Centroid and Peak 

Chandra and HST images of IDCS J1426.5+3508 are 
shown in Figure The first column shows Chandra 
images and contours for Gaussian-smoothed (FWHM = 
3", upper panel) and adaptively smoothed {lower panel) 
0.5-2.0 keV X-ray images, respectively. The 2nd column 
shows these contours overlaid on color optical/IR HST 
images, a full description of which are given in (Mo et al. 
in prep.). The upper right panel is zoomed in to the 
cluster core to better show the brightest cluster galaxy 
(BCG) and the giant gravitational arc. 

We use the centroid measured within a 250-500 kpc 
annulus as the cluster center, which tends to provide less 
biased estimates of the global properties for unrelaxed 
clusters. This choice of center, (ax, Sx) = (14:26:32.6, 
+35:08:25), is within 4" of the X-ray peak , within 5'' 


(5.3 a) decrement 
IDCS JI426.5+3508 has a giant gravitational arc, from 
which a mi nimum mass of Msnn a rr ^ 2.0 x 10 ^"^ Mq was 
estimated (Gonzalez et al. 2012). Although these inde- 


of the BGG position ( Stanford et al.||2012 ) and within 
28" of the SZ centroid (jBroHwi^W^rl 3012). Given the 
uncertainty in the SZ centroid (« 35'"), this positional 


Brodwin et al. 2012 1 . Finally, offset is not statistically significant, 
las a giant gravitational arc, from 


pendent mass measurements are all in good agreement, 
indicating that IDCS J1426.5+3508 is a very massive, 
relaxed cluster at z = 1.75, the Lx-based X-ray mass 
was based on only 53 counts and is highly uncertain. 

In this paper we present deep new Chandra obser¬ 
vations from which we measure the ICM properties of 
IDCS J1426.5+3508. In (j^we describe the data used in 
this analysis. In )|^we present the flux, energy spectrum 
and gas density profile of IDCS J1426.5+3508, along 
with the quantities we derive from these direct measure¬ 
ments, including the luminosity, temperature, gas mass 
and metallicity. Using standard scaling relations from 
the literature, we estimate M 500 for IDCS J1426.5+3508 
from four different X-ray estimators in ® and compare 
these to complementary SZ- and lensingTbased mass es¬ 
timates. We also compute gas fractions for each of these 
halo mass estimates and compare these with the value 
predicted from low-redshift clusters. In ^ we place 
IDCS J1426.5+3508 in the context of the known z > 1.5 
galaxy cluster population and discuss the evolutionary 
state of its ICM. Finally, we present our conclusions in 
We use rim = 0.27, Da = 0.73 and Hq = 70 km s“^ 
Mpc“^ throughout. 

2. DATA 

The X-ray data used in this work were acquired via 
Chandra proposal 148000534 (PI: Brodwin). A total ex¬ 
posure of looks was acquired over two pointings (OB- 
SIDs 15168 and 16321). This exposure time was chosen 
to obtain 500 X-ray source counts, but due to the evolv- 


3.2. Point Sources 

The incidence of AGN in clusters has been shown 


2015 


(301 


to increase rapidly with redshift (Eastman et al. 2007 
Gal ametz et al.|[2009[|Martini et al.||2013 
This is natural consequence of the 


Alberts et al. 

consequence ot tlie Brodwin et al. 
merger model discussed above, where the mergers 
that induce starbusts also fuel the AGN that provide the 
eventual quenching. 

Point sources were identified and masked using an 
automated routine follo wing the wavelet decom position 
technique described in Vikhlinin et al. (1998). There 
are two bright X-ray poi nt source near the clust er core, 
previously discussed in Stanford et al. (2012). The 
northern one is a QSO in the cluster, with a flux of 

erg s“ 


'S'n 5-2 ~ 1.85^n,7 X 10 erg s ^ cm The south- 


0 . 5-2 ~ r.ou_p 37 X iU erg s " cm 
ern one is a bright ra dio source liste d in both the NVSS 
([Condon et al. 1998[) and FIRST ([Becker et al.[[1995l 
catalogs. It has a tlux of ~95 mJy at 1.4 GHz, and was 
found to have a 31 GHz flux of 5.3 ± 0.3 mJy in the SZ 
analysis of Brodwin et al. (20121. Here we measure a flux 
of S' 0 . 5-2 « 5.68l(Q;gg X 1(1“^^ erg s“^ cm“^. These point 
sources, along with several others at larger clustercentric 
radii, are masked and do not affect this analysis. 


3.3. Counts, Flux and Luminosity 

We measure 401 net counts (0.5-6 keV) from 
IDCS J1426.5+3508, after point source masking and 
background subtraction. The flux in the soft band is 


S' 0 . 5-2 = (2.2 + 0.6) X 10-14 erg s"! 


cm corresponding 


ing effective area at low energies (Marshall et al. 2004 
O’Dell et al.|2013' ) we obtained slightly fewer than antici¬ 


pated. The data, obtained with ACIS-I, were cleaned for 

® refers to the mass within r^, the radius at which the 

mean overdensity is A times the critical density, as inferred via 
mass proxy A. 


to a luminosity of + 0 . 5-2 = (3.6 ± 0.5) x lO^^ erg s”!. 
The bolometric luminosity over the 0.01-100 keV energy 
range is m easured from the best-fit model described be¬ 
low ((3.4) to be Lx,boi = (12.8 ± 1.1) x lO^^ erg s”!. 
These quantities are all measured within rsoo.Xx = 530 + 
10 kpc, as determined from the co re-excised Y x mea- 
surement described below ((|3.6[). [Pratt et al.[ ([20091 
showed the bolometric luminosity to be a lower scatter 
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100 kpc I 100 kpc 



Fig. 1.— Upper left: Gaussian-smoothed (FWHM = 3^^) C/iandra 0.5-2.0 keV image of IDCS J1426.5+3508. Contours correspond to 0.02, 
0.05, 0.1, 0.18, 0.22, 0.24, 0.25 and 0.45 counts per 0.492^^ xO.492^^ pixel. The location of the core as traced by the ICM peak is indicated. 
Upper right Pseudocolor ACS/F606W, ACS/F814W and WFC3/F160W HST image of the central region of IDCS J14 26.5+3508, with 
the X -ray contours from the upper left panel overlaid. The giant arc, boxed, is visible ~125 kpc north of the BCG (see |Gonzalez et al.| 
[20^. Lower left: Adaptively smoothed image of Chandra 0.5-2.0 keV image of IDCS J1426.5+3508, showing the large-scale structure 
trom which the cluster centroid was measured. The contours correspond to 0.011, 0.016, 0.022, 0.038, 0.065, 0.085 and 0.15 counts per 
0.492" X 0.492" pixel. The peak emission in the core is offset by ~30 kpc from the cluster centroid. Two bright point sources discussed 
in the text are indicated. Lower right: The HST image with the same scale and contours as the adaptively smoothed X-ray image. 


mass proxy than Lq 5 _ 2 . Using their scaling relation, we 
find Msoo.ix.boi = (2.8 + 0.3) x Mq. 

3.4. Temperature 

The 0.5-6. 0 keV ACIS-I spect ra within r 25 oo 


(^0.45r5oo,Yxi Vikhlimn et al. 2006) for each OBSID 
are shown in 1* igure TEIs aperfure was chosen to 
maximize signal to noise, although we also consider a 
core-excised (0.15-l)r5oo annulus below. The spectro¬ 
scopic temperature was measured by modeling the X-ray 
spectrum using a combination of an absorbed, optically 
thin plasma (phabs x APEC), an absorbed hard back¬ 
ground component (phabs X bremsS; kT = 40 keV), 
and a soft. Galactic background component (apec; kT 
= 0.15 keV). Foreground and background models were 
constrained by fitting simultaneously to an off-source re¬ 


gion within the same field of view and to the on-source 
region. The Galactic hydrogen column density, Nh, was 
set to the wei ghted average from t he Leiden-Argentine- 
Bonn survey (iKalberl a et al.||2005|). The s ource redshift 
was fixed to z = 1.75 ( btantord et al.|2012 ). The reduced 
of this fit is 0.95, demonstrating that the data are very 
well described by the simple, single-temperature plasma 
model. While the goodness-of-fit is excellent, the indi¬ 
vidual parameters (e.g., kT, Z) are not well constrained, 
as discussed below. 

The average temperature within r 25 oo is kT 25 oo = 
6.2tf 0 keV. This global (not core-excised) temperature, 
while not optimal as a mass observable due to inclusion 
of the cluster core, is useful for comparison to lower res¬ 
olution measurements from other facilities or low SNR 
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Chandra observations for which all the counts are re¬ 
quired to measure the temperature. For completeness 
we also report a core temperature, within r < 0.15 rsoo, 
of fcTcore = 7.31:11 keV. 

We measured the core-excised spectroscopic temper¬ 
ature over (0 15 — 1) r^nn using the Tjf-Msoo scaling 
relation from |Vikhlinin et al. (2009a) to estimate M 500 , 
iteratively updating the radius (rsooj until it converged. 
We used t he pipeline descri b ed in Benson et al. (in 
prep.) and McDonald et al. ( 2013), which closely fol - 
lows the procedures described in Andersson et al. (2011). 
We find kT^oo = 7.6^®;g keV, in good agreement with 
the global value given above. This temperature corre¬ 
sponds to a mass of M^qq^Tx = ^ 

^ 500 , = 560 ± 20 kpc, where the subscript refers to the 

scaling relation from which the physical quantity was de¬ 
rived. The error on the core-excised is larger than 

that on kT 25 oo as roughly a third of the signal is removed 
with the core, and the noise is significantly increased 
by moving from r 25 oo to rsoo- Nevertheless, tempera¬ 
ture measurements and confidence level determinations 
remain unbiased down to count levels well below those 
in this work (e.g., Churazov et al.|[T9M|. 
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Fig. 2.— Upper panel. Energy spectrum of IDCS J1426.5-E3508 
interior to r2500. The gray points are the two separate observations 
to which the models are fit, and the black points show the stacked 
spectrum. The global temperature in this central region, based 
on the overplotted best-fit model, is kT2500 = 6.2^)'® keV. The 
expected location of the Fe Ka line is indicated, though we do 
not detect it. The 1 a limit on the central metallicity is Z2500 < 
0.33 Zq, as dicussed in j ]3.7| Lower panel. Residuals of the best-fit 
model. 


3.5. Gas Mass 

Similarly, Mn. 5 nn was derived via t he /g-Msoo scaling 
relation from Vikhlinin et al. (2009a). The derivation of 
the gas den sity prohle, described in detail in McDonald] 


et al. (2013), involves measuring the X-ray surface bright- 
ness in the rest-frame energy 0.7-2.0 keV as a function 
of radius. The resulting surface brightness profile, shown 


in Figure [3 
beta moc 


{upper yan el ), is fit with a project ed double- 

(2006). In con- 


following jVikhlinin et al. 
verting from electron density to gas density, we assume 
Pg = mpUeAlZ, where A = 1.397 is the average nuclear 
charge and Z = 1.199 is the average nuclear mass. We 
integrate the deprojected gas density profile within an 
initial radius of estimate the total enclosed gas 

mass, and then use this to derive M 500 via the /g-Msoo 
relation. This process is iterated until convergence, lead¬ 
ing to a revised estimate of rsoo, Mg = 500 ± 10 kpc and 


M, 


500, Mj 


= 2.3 


- 1 - 0.7 

- 0.5 


X 10^^ Mq, consistent with values in¬ 


ferred from the Tx-M scaling relation. The gas mass is 

Mg,500 = 2.5t° ® X 1013 Mtg. 


''- 0.6 


E 





Fig. 3.— Upper panel: X-ray surface brightness profile for 
IDCS J1426.5-^3508, derived using the large-scale X-ray centroid 
as the center. The gray curve represents the best-fit projected 
double-beta model. Lower panel: X-ray surface brightness profile, 
centered on the X-ray peak. The gray curve is the same as in the 
upper panel. There is evidence of a dense core on scales of ~ 20 
kpc, which appears slightly offset from the large-scale centroid (see 

Fig.g. 


3.6. Fx 

The product of the core-excised temperature and gas 
mass, referred to as Yx, approximates the total thermal 
energy in the cluster. It has been shown in simulations 
to be a low-scatter mass proxy t hat is independent o f 
the dynamical state of the cluster ( jKravtsov et al.|2006 ). 
Following the same approach as tor k'l '500 and Mg, 500 , 
we estimate Yx, 5 oo iteratively, adjusting rsoo to satisfy 
the Yx-Msoo relation. We find rsoo.Xx = 530 ± 10 kpc 
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and Msoo.Yx = ^ Within rsoo.Vx; 

measure I500 = l-S^o'e ^ 10^"^ ^^0 keV. 


TABLE 1 

Properties oe IDCS J1426. 5+3508 


3.7. Metallicity 

Interior to r 25 oo and r^oo (non core-excised), we 
measure Itr (2(t) upper limits on the metallicity of Z 2500 
+ 0.33 Zq (0.47 Z 0 ) and Z 500 + 0.10 Zq (0.18 Zq), 
respectively. As Figure [^illustrates, there is no obvious 
detection of the Fe Ka emission line, such that, within 
^ 2500 ) these data are consistent with anywhere from 
“typical” to a complete absence of metals. Over larger 
radii we may be observing a marked lack of metals, with 
the typical average metallicit y of low-redshift clusters 
being Z 500 ~ 0.3 ± 0.1 Zq (De Grandi and Molendi 


clusters (Rosati et al. 2009 

Santos et al. 

2012 

De 

Grandi et ai. 

2U14). However, deeper data are needed 


to properly constrain the metallicity in this system, and 
determine whether it is, indeed, metal-poor. 


4. TOTAL MASSES FOR IDCS J1426.5+3508 

4.1. M 500 from X-ray Scaling Relations 

Using the measured luminosity, temperature, gas mass 
and + 500 , we reported in the previous section several 
complementary estimates of the total mass within rsoo. 
These masses, along with a host of physical parameters, 
are listed in Tablejl] Had we assumed no evolution rather 
than a self-similar evolution, our Fx-based mass would 
be ^ 50% higher. 

These masses all show a high level of consistency, indi¬ 
cating IDCS J1426.5+3508 is relaxed, and offering good 
support for the hydrostatic approximation that is under¬ 
lying all these scaling relations. A more stringent test is 
provided by comparisons with complementary, indepen¬ 
dent mass probes. 


4.2. M 500 from Complementary Measurements 

IDCS J1426.5+3508 is unique in that it is the high¬ 
est redshift cluster, by far, for which independent mass 
measurements are available in the X-ray, SZ and from 
strong gravitational lensing analyses. This rare conflu¬ 
ence is only possible due to its extreme mass — it is the 
most massive cluster known at z > 1.5 — and the fortu¬ 
itous (and surprisi ng) presence of a ba ckground galaxy 
for strong lensing (|Gonzalez et al.||2012|). 

Brodwin et al.j J 2012 p reported a strong (> 5cr) SZ 
detection at 31 GHz with CARMA. The resulting mass 
measurement, Msoo.sz = (2.6 ± 0.7) x 10^^ Mq, is in 
excellent agreement with all of the X-ray mass measures 
prese nted above. Relative to the Yg^-M scalin g relation 
from lAnder sson et al. ([2011 1 that was used in Brodwin 
et al.fl 2012 ), the normalization is expected to increase by 
~25% in Ysz, corresponding to a ~15% increase in mass 
(Benson et ah, in prep.). The increase is primarily due 
to two factors — a statistical shift resulting from a much 
(^5x) larger SZ and X-ray cluster sample, and updated 
mass no rmalization from more recent weak-lensing obser¬ 
vations (jHoekstra et al. [2015). The close agreement of 
the SZ rnass tor IDCS J 1426.5+3508 with the X-ray mass 
measures will not be aff ected by this modest change. 


Gonzalez et al. (2012^ reported the discovery of a giant 
gravitational arc, visible in Figure about 15" N of the 


Property 

Value 

Unit 

z 

1.75 


X-Ray Centroid 

(14:26:32.6, +35:08:25) 


X-Ray Peak Position 

(14:26:32.9, +35:08:26) 


SZ Centroid 

(14:26:34.0, +35:08:03) 


BCG Position 

(14:26:32.95, +35:08:23.6) 


^500, 

560 ± 20 

kpc 

rsoo, Mg 

500 ± 10 

kpc 

r500,Yx 

530 ± 10 

kpc 

Counts 

420 


Count Rate 

4.14 xlO-3 

s-l 

So. 5-2 

(2.2 ±0.6) X lO-i-^ 

erg s~^ cm“^ 

So. 5-2 

(3.6 ± 0.5) X 10*^ 

erg s“^ 

Sx. bol 

(12.8 + 1.1) X lO-^^ 

erg s“^ 

kTcmC 

7 0+2.7 

'•'^-1.3 

keV 

kT2500 

6 2^^'^ 

'^•^-1.0 

keV 

kTsoo 

7 

'•'^-1.9 

keV 

Mg,500 

2.5t° g X 10^3 

Mq 

Y 500 

l-9to 6 X 1014 

Mq keV 

■^2500 

< 0.33 (lo-) 

Z0 

•^500 

< 0.10 (lo-) 

Z0 

M 500 , Yx 

2.6/J;| X 1014 

Mq 

^500,Tx 

3.3ti-2 X 1014 

Mq 

M 5 OO, Mg 

2.3t)J/ X 1014 

Mq 

M 500 , Lx, bol 

(2.8 ±0.3) X 10i4 

Mq 

^500, SZ 

(2.6 ±0.7) X 10i4 

Mq 

^500, arc^ 

(1.9/°/ - 2.6/°;®) X 1014 

Mq 


^ The core temperature is measured within r < 0.15 rsoo. 

^ This mass ranges account for possible source redshifts beween 4.5 < 
2 < 6 , as described in the text. 


BCG at a clustercentric radius of ^125 kpc. Attempts to 
secure a spectroscopic redshift for the source galax y were 
unsuccessful, thou gh from the available photometry |Gon-] 
zalez et (20121 constrained the source redshift to the 


range 2 Zs 6 . With a non-detection in subsequent 
deep (AB ^ 28, 10 a) HST/ACS F606W imaging, new 
data described in (Mo et al. in prep.), the redshift con¬ 
straint is now refined to 4.5 ^ z ^ 6. The strong lensing 
mass directly measured within the arc radius ranges from 
(6.9 ±0.3) X 1013 Mq for z,, = 6 to (8.5 ±0.3 ) x 10^3 Mq 
for z,^ = 4.5. Extrapolating to rsoo using the [Duffy et al.| 
(20081 mass-concentration relation, this corresponds to 
a mass range between M500, arc = l-S^o's x Mq and 
Afsoo,arc = X lOi^ Mq, in good agreement with 

all the ICM-based masses. 

A weak lensing analysis of IDCS J1426.5±3508 is un¬ 
derway (Mo et al. in prep.). A shear signal is detected in 
Cycle 20 HST images, consistent with expectations for a 
cluster this massive, even at z = 1.75. 

Figure [^ {upper panel) compares all of the mass mea¬ 
sures described in this paper, and shows the uniformly 
excellent agreement among them. This confirms that, 
despite its extreme redshift, IDCS J1426.5±3508 is by 
all measures a relatively evolved, relaxed cluster. We 
take the low-scatter Yx-based mass as our best esti¬ 
mate of the halo mass of IDCS J1426.5±3508, Msoo.Vx 
= X 10 ^"* Mq. 


4.3. /gas 

As massive galaxy clusters assemble, they are expected 
to retain baryon fractions {fb) close to (but slightly 
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Fig. 4.— Uvv er yanel: Masses for IDCS J1426.5+3508 from the X-ray measu rements described herein (filled circles), along with com¬ 
plementary SZ ||Brodwin et al.||2012l) and strong leasing (iGonzalez et al.||2012|| masses (filled squares). For the latter we plot the two 
masses corresponding to tne extreme source redshift values, as described m tne text. We take t he Yx-based mass our best estimate (green 
horizontal line and error range), as it has the lowest intrinsic scatter ||Kravtsov et al.||2006|l. The remarkable agreement among these 
independent methods demonstrates that IDCS J1426.5+3508 is a remarlbibly mature, evolved cluster even at z = 1.75. Lower panel: fgas 
measurements for IDCS J1426.5+3508 using each of these halo mass measurements, fgas measurements with independent M^^soo and halo 
mass measurements have filled symbols; those affected by covariances between the gas and total mass are shown with op en symbols. The 
expect ed fgas value for IDCS J1426.5+3508, shown as the red line with 1 a errors, is taken from the fqas-M^QQ relation of|Andersson et al.| 
( |2011[ |. The observed value of fgas is in good agreement with the predicted value for all mass probes. 


lower th an) the universal value, with r elatively low scat- 
ter (e.g., Kravtsov and Borgani 2012 1 . Observationally, 
the gas rnass traction [fgas) is highest in the most mas¬ 
sive clusters, with a weakly decreasing fraction to lower 


mass es (e.g., Vikhlinin et al. 

rnll. 


2009a Andersson et al. 


We calculate fgas using each of the masses derived 
above, integrating the gas mass out to the rsoo value ap¬ 
propriate to each mass proxy. These are plotted in Fig- 
ure Id] {lower panel). The fgas measurements for which 


This provides additional evidence of the maturity of 
IDCS J1426.5-1-3508 and bolsters our confidence in our 
use of scaling relations at this extreme redshift to esti¬ 
mate its mass. The recent measurement of the gas frac¬ 
tion in XDCP J0044.0-2033 at z = 1.58 (Tozzi et al. 
2015|), fgas= 0.08 + / — 0.02, is slightly lower than mean 
ot the present measurements, though in agreement within 

1 (T. 

5. DISCUSSION 

the Wlg_ 50 o and halo mass measurements are indepen- 5.1. Comparison With Other Distant, Massive Clusters 


dent are shown as filled symbols, while those affected by 
covariances between the gas and total mass have open 
symbols. We compare the measured fgas valu es with the 
value predic ted from the redshift-independent |Andersson| 
et al. (2011 1 relation. This is shown as the solid red line, 
where the error r ange includes the e rrors from fgas-M^Qo 
fit parameters in Andersson et al. (2011), as well as the 
mass error in M 5 oo(ly 7 r"^^IiicIi~we took as the value of 
M 500 for this relation. 

The observed value of fgas from every available mass 
proxy is completely consistent with the expected value 
for a cluster with the mass of IDCS JI426.5-1-3508. 


IDCS J1426.5-1-3508 was the first relaxed, massive 
galaxy clus ter to be confirmed at z > 1.5 (|Stanford 


et al. 2012). The present analysis of its X-ray prop- 
erties conhrms t he SZ and lensing masses r eporte d in 
Brodwin et al. (2012) and Gonzalez et al. (2012), re- 
spectively. 'i'he tornier paper also demonstrated that 
IDCS J1426.5-1-3508 is an evolutionary precursor to the 
most massive known clusters at all redshifts. Sev¬ 
eral new z > 1.5 clusters have subsequently been re¬ 
ported, but none are as massive, and hence rare, as 
IDCS J1426.5-1-3508. Although the probability of detect¬ 
ing a cluster this massive in the 9 deg^ IDCS area is very 
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small (< 1%), Brodwin et al. (2012) demonstrated that 
its existence poses no threat to ACDM. That paper also 
estimated that the SPT should expect ^ 2.4 such clus¬ 
ters over their entire 2500 deg^ s urvey. In the final SPT 
catalog paper, |Bleem et al.| (|2015|) reported three clusters 
with strong decrements that do not yet have optical/lR 
confirmations. Though not confirmed to date, current 
photometric limits suggest that these three clusters lie 
at z = 1.7 ±0.2. 



low the identical mass measure for IDCS J1426.5±3508, 
^500, Yx= X 10^"^ Mq, as well as all the other 

ICM-based (i.e. X-ray and SZ) mass measures in Table 

m 

A direct comparison of tem peratures is not straight¬ 
forward, as |Tozzi et al.| (|2015|) do not quote a value for 
72500 or Tsoo^ 'I'hey instead measure a spectroscopic 
temperature of kT = 6.7lJ;9 keV at a radius of 375 
kpc. As the raoo values for both clusters are quite sim¬ 
ilar, we measure the temperature at this metric radius 
in IDCS J1426.5±3508 in order to make a meaningful 
comparison. We find fcTayakpc = 7.3^2’g keV, confirming 
that IDCS J1426.5±3508 is likely the hotter and more 
massive cluster. 


Newman et al. (2014) presented the spectroscopic con¬ 


firmation of JKOS 1)41 at z = 1.80, for which Andreon 


et al. (2014) report a mass, based on its X-ray and optical 
properties, in the range M500 ~ (1 —2) x lO^'^ Mf^. Given 


the n on-detection in deep SZ imaging (Culverhouse et al. 


20101, the mass is likely at or below the lower end of that 


range and thus less massive than IDCS J1426.5±3508. 
Finally, cluster c a ndida t es reported at 1 . 6 5, ^ < 2 
■ Go- 


by Papovi c h et a l. (|201 0), Z eimann et al.| ( 2012|) , 
bat et al. 1 2011 ” 20131 ana iVlantz et al.| (2014) 


are 


all considerably les s massive than IDCS J 1426.5-1-3508. 


Mantz et al. (|2014|) describe the most massive of these, 

' l-2x 1014 Mq, 
1.9. Despite this 


XLSSU J021744.1-034536 with M500 
at a photometric redshift of Zphot 
redshift uncertainty, which strongly affects mass proxies 
in the X-ray, we can say with certainty that this clus¬ 
ter has a lower mass than IDCS J1426.5±3508. This 
assertion is based on a comparison of the spherically 
averaged dimensionless Comptonization, Psph.soO) mea¬ 
sured with CARMA for both clusters. The SZ ma,ss 


proxy M PC (Msph.aoo 77yA(z)2/3)3/5 (e.g. 


^ spn,ouuV^/ / 

et al.|[2012|, where Da is the angular diameter 


Marrone 


distance 

and D{z} IS evolution of the Hubble parameter, is weakly 
dependent on redshift. With Psph.soo being a factor of 2.6 
higher in IDCS J1426.5±3508, we find its mass is larger 
by a factor of ^1.9. 


5.2. Dynamical and Cooling State of the ICM 

In Figure we show the projected surface brightness 
profile for two different choices of center: the large- 
scale (250-500 kpc) centroid of the X-ray emission and 
the X-ray peak. Using the former definition, we fit a 
projected beta model to the data, finding no evidence 
for a central surface brightness excess. However, we 
show in the lower panel of Figure that the peak of 


the X-ray emission — which lies ^30 kpc from the cen- 
troi d — represents a sign ificant overdensity. Follow¬ 


ing Vikhlinin et al.| (2007), we measure the cuspiness 
of the peak-centered surface brightness profile, finding 
a = (dlogpg/dlogr) |o.041-500= 0-^2 ± 0. 09. Such a high 
cuspin ess at high redshift is rare. Indeed, |Vikhlinin et al.| 
(2007) found no such systems at z > 0.5 in their 400 deg^ 
survey. Two similar systems, albeit at much lower red¬ 
shift (z ^ 1. 1), have been identified — one by McDonald 
^_^y_|2013) i n a sa mple of 83 SPT clusters, and one by 


Santos et al.| (20121 in the WARPS survey. 

I'he offset of ~30 kpc between the X-ray centroid and 
the dense core suggests that this cluster has undergone 
a recent interaction, and that the cluster core is sloshing 
about the potential minimum. Such an offset ou ght to 


remain visible for < 500 Myr after any interacti on (Asca- 
sibar and Markevitch 2006 ZuHone et al.|20f0). Consid¬ 


ering that this system is being observed when the Uni¬ 
verse was only 3.8 Gyr old, and that it had to assemble 
rapidly to achieve such a high mass at such early times, it 
is unsurprising that it retains an imprint of this hurried 

growth in it s core. _ 

Following jMcDonald et al. (2013), we calculate a 
pseudo-deprojected entropy, using the deprojected gas 
density profile (§3.5), an aperture temperature, and as¬ 
suming the X-ray peak as the center. We find Kq ~ 20 
keV cm^, corresponding to a cooling time of ~ 160 Myr. 


These properties are typical of cool co re clusters (Hud¬ 


son et al. 2010 McDonald et al. 2013), suggesting that 


such systems can form very early in the cluster lifetime. 

The fact that the dense core appears elongated (Figure 
suggests that, rather than a traditional cool core, we 
may be observing a dense infalling group. This is consis¬ 
tent with many of the observed qualities, including the 
low entropy, the offset from the potential minimum and 
the non-symmetric morphology. Conversely, the near! 
coincident BCG and X-ray peak (within ^ 3"; Figure [I 
suggests that this may be an offset core rather than an 
infalling group. Distinguishing between these two sce¬ 
narios requires a deeper X-ray follow-up observation. 


6. CONCLUSIONS 

We have presented a deep 100 ks Chandra observation 
of IDCS J1426.5±3508 at z = 1.75, the most massive 
cluster discovered at z > 1.5 from any method. We mea¬ 
sured the luminosity, temperature and gas mass, from 
which we derived halo mass estimates from the Tx-M, 
/g-M, Yx-M, and Lx-M scaling relations. These all 
show excellent consistency and are in remarkable agree¬ 
ment with independent SZ and strong lensing masses. 
Similarly, the gas mass fractions for all these mass prox¬ 
ies were found to be in good agreement with each other 
and with the value predicted from low redshift clusters. 

The bolometric luminosity is Ax, boi = (12.8 ± 1.1) x 
1944 gj.g g-i^ from which we estimate a mass of 
M500,Lx.boi = (2.8 ± 0.3) X 10^4 Mq. We measure a 

central temperature within r25oo of kT 2 ^m = 6.2j'ji[ Q keV 
and a core-excised temperature within rsoo, = 560±20 
kpc, used in mass scaling relations, of fcTsoo = 7.6^1(9 
keV. This results in a mass of M500, Tx = 3-3ii;2 ^ 1^44 
Mq. We find no evidence for metals in the ICM of this 
system, placing a 2a upper limit of ■^500 < 0.18^0, sug¬ 
gesting that this system may still be in the process of 
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enriching its intracluster medium. 

We measure a gas mass of Mg^soo = 2.5lg;g x 10^^ Mq, 
from which we infer an Mg-based halo mass of M 500 , Mg = 
2 . 3 to ;5 ^ -^ 0 - From the gas mass and core-excised 

temperature, we find Y 500 = ^ -^0 FeV, 

from which we derive our lowest-scatter estimate of the 
mass of IDCS J1426.5-h3508, M 5 oo,Yx= x 10^"* 

Mq. All the X-ray masses for IDCS J1426.5-1-3508 are 
in agreement, as are the SZ- and lensing-based mass from 
previous analyses. 

The cluster has a dense, low-entropy core, offset by 
~30 kpc from the X-ray centroid, which makes it one 
of the few “cool core” clusters discovered at 2 : > 1, and 
the first known cool core cluster at 2 : > 1.2. The off¬ 
set of this core from the large-scale centroid suggests 
that this cluster has had a relatively recent (< 500 Myr) 
merger/interaction with another massive system. We 
measure a central entropy of Kq ~ 20 keV cm^, indi¬ 
cating the cool core may have a very rapid cooling time 
of ^ 160 Myr. 

IDCS J1426.5-1-3508 is the first cluster at z > 1.5 to 
have all these independent mass measurements. In ad¬ 
dition to being the most massive known cluster in this 
redshift regime, IDCS J1426.5-1-3508 also has a remark¬ 
ably relaxed ICM suggesting a very early and rapid for¬ 
mation. Deeper follow-up X-ray observations are essen¬ 
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